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Abstract 

We have developed a general approach for the calculation of the single molecule polarization 
correlation function C(t), which delivers a correlation of the emission dichroisms at time and 
t. The approach is model independent and valid for general asymmetric top molecules. The 
key dynamic quantities of our analysis are the even-rank orientational correlation functions, the 
weighted sum of which yields C (t) . We have demonstrated that the use of non-orthogonal schemes 
for the detection of the single molecule polarization responses makes it possible to manipulate the 
weighting coefficients in the expansion of C(t). Thus valuable information about the orientational 
correlation functions of the rank higher than second can be extracted from C(t). 



I. INTRODUCTION 



The formulation of the polarization-sensitive spectroscopy in terms of the orientational 
correlation functions (OCFs) has opened up a possibility of the unified description of various 
signals and clarified information content of different spectroscopic techniques." As has been 
demonstrated, all large variety of different polarization signals can be described by the OCFs 
of the first and second rank.— More specifically, one can "measure" either OCFs in the 
time domain, or their time derivatives, or their Fourier spectra, or their integral relaxation 
times. It is this unification which has made it possible to compare the results of different 
measurements and to learn about the mechanisms of molecular reorientation, both in the 
gas phase and in the condensed phase. The contemporary nonlinear (third order) ultrafast 
polarization spectroscopy is interpreted in terms of the three-time correlation functions 
(CFs) of the dipole moments or polarizability tensors involved.- However, in practices, their 
influence on a signal is normally accounted for by a static averaging, or the three-time 
polarization CFs reduce to the standard OCFs of the second rank, due to either the strong 
optical dephasing or the shortness of the laser pulse on the rotational dynamic timescale. 

All the written above pertains to spectroscopy of ensembles, in which the measured 
response is averaged over many single-particle contributions. The situation with the single 
molecule (SM) spectroscopy in general,- 1 ^*^ and the SM polarization spectroscopy in 
particular,-*^ is very different. The SM signal delivers a response of an individual system, 
whose time dependence reports about fluctuations caused by the system " nanoenvironment" . 

Ideally, one would prefer to measure the three-dimensional orientation of the emission 
dipole moment (s) in real time. Indeed, there exist several schemes which allow us to do 
that^^i^i^Lillill These techniques require, typically, many photons to get a good signal- 
to-noise ratio. This obstacle restricts the length and the time resolution of the recorded 
signal. Fourkas has suggested that the time evolution of the emission dipole can be deter- 
mined "on the fly", by detecting SM emission along three different polarization directions.— 
Hohlbein and Hiibner have recently implemented this method.— 

So far, many SM experiments have been designed to detect the in-plane projection of the 
SM emission along two different (mutually perpendicular) polarization directions. Thereby 
long enough transients can be measured, which permit of the reliable calculation of SM 
CFs. The SM dichroism CF is the key dynamic quantity which is delivered by the polariza- 
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tion SM fluorescence microscopy.— i22j2L22a22a2i This CF is much more complicated than its 
counterparts which describe polarization responses in ensemble measurements. 

The SM dichroism CF is the core object of the present study, which has been inspired 
by the recent papers.— 1 ^ Our aim is threefold. First, we generalize the approach, which has 
been developed in22>2i f or spherical top molecules within the small-angle rotational diffusion 
model, to asymmetric top molecules and beyond any particular model of molecular reori- 
entation. The key dynamic quantities here are the even-rank OCFs, which can either be 
evaluated within any model of the molecular reorientation available in the literature,- 1 *^ or 
simulated on a computer.— Second, we discuss a possibility of gaining additional knowledge 
about molecular reorientation by utilizing non-orthogonal schemes for the detection of SM 
polarization responses. Third, we demonstrate that a valuable information about the OCFs 
of the rank higher than second can be extracted from the SM CFs. Until recently, such an 
information was available only through the computer simulations and model calculations. 

II. SINGLE MOLECULE DICHROISM 

The key quantity in the SM polarization microscopy is the dichroism 



Here Ii are the intensities of the light emission which are detected at two different (usually, 
mutually perpendicular) polarizations e*, i = 1, 2. Within the oscillator model, the excita- 
tion and emission processes are independent, so that the emission intensity is proportional 
to the product of two probabilities,- 1 ^ 



and the SM dichroism is independent of the absorption cross-section a a . The absorption 
probability is given by the square of the scalar product of the polarization of the absorbed 
light e a and the absorption dipole moment /i a : 28 



Va(0)ai{t), 



(2) 




(3) 



The emission cross-section is given by a similar expression,— 




(4) 



3 



in which the polarization of the emitted light, e, is explicitly given by the equation 

£ = kx/jxk = /i - k(/ik), (5) 

k being the unit vector along the propagation of the light beam, and \i being the emission 
dipole moment.™ 

When the fluorescence is collected from an ensemble of molecules, molecular contributions 
add up incoherently. Thus, in order to get the total emission intensity, one has to average 
Eq-dSD over all possible orientations of the wave vector k. This yields (e) k ~ fi, so that the 
averaged emission cross-section becomes similar to its absorption counterpart: 

<n ~ (e,/i) 2 . (6) 

Eq.(j2j), in conjunction with Eqs.flS]) and ([H]), embody the standard starting point for the 
calculation of the intensity of the polarized emission in ensemble measurements, (Ii(t)) ~ 
(a a (0)ai(t)), see Refs.^ 

If we collect emission from a single molecule, the above incoherent averaging procedure 
is not legitimate any longer. Since the wave vector k is specific to a photon which has been 
emitted by the molecule but the direction of k is unknown, we have to average the emission 
probability (jl]) over all those k which can be collected by the objective. Let s be the unit 
vector along the axis of the objective (s_Lej) and i?k, 0k be the spherical angles of the unit 
vector k, so that ks = cos($k)- We also introduce the light-collection cone angle 

X = arcsin(iWi/n), (7) 

NA being the numerical aperture of the objective and n being the refraction index of the 
medium in which the sample is embedded (Fig. 1 clarifies the meaning of the introduced 
quantities). If we further assume that the emitting molecule is located at the focal point 
of the ideal and polarization-preserving objective which collects the emission, then the k- 
averaging of the emission probability fll]) yields 

,2tt fX 



/ dcj) I wv&d&ai ~ p + g(s/i) 2 + 2(e^) 2 . (8) 

Jo J-x 



The numerical parameters p and q are uniquely determined by the collection angle 

p—(A + B)/C—l, q=l-B/C, 



v . 11,18,23,24,32 



where the quantities A, B, C are defined in the standard way: 

A = 8 - 12 cos(x) + 4 cos 3 (%) , (9) 

B = 6cos(x) - 6cos 3 (x), (10) 
C = 7 - 3cos( X ) -3cos 2 ( X ) -cos 3 ( X ). (11) 

Eq. (jSj) generalizes slightly its counterparts from,— il&i2&2£i22 allowing for an arbitrary direction 
of the polarizer ej_Ls in the laboratory frame. 

Fig.2 shows the behavior of the coefficients of p(x) and q(x)- Both of them increase 
monotonically with \: p(0)= q(0)= 0, while p(ir/2) = 1/7 and q{n/2) = 1. For every Xi 
q(x) > p(x)- Eq.fJHJ) reduces to Eq.(jSJ) in the limit of small collection angle (x — ► 0). In 
this case, the beams with k||s are collected only, the k-dependent portion of Eq. (J5J) does not 
contribute into Eq.(jl]), so that Eq.flS]) reduces to (JOJ) since both p — > and q — > 0. 

After the insertion of Eq.® into Eq.flT]), one gets the general expression for the SM 
dichroism: 

(e!/i) 2 — (e 2 //) 2 

X = p + g(s /U )2 + (e 1/U )2 + (e 2/ i)2- (12) 
If NA = 0, then p = q = and the dichroism attains the familiar form 

v _ (ei/i) 2 - (e 2 /i) 2 
A -(erf + (e 2 ^- 

If we employ the orthogonal signal detection scheme (ei_Le 2 ), then Eq. (fT2|) simplifies to 

l+p- (1 -g)(s/i) 2 v ; 

If, additionally, A^A is high (\ ~* 7r /2, p(ir/2) = 1/7 and q{ir/2) = 1), then the denominator 
in Eq. f|T4l becomes constant and 

X = I{(e 1/U ) 2 -(e 2/ u) 2 }. (15) 

III. ORIENTATION AL CORRELATION FUNCTIONS 

Let O denote collectively the set of three Euler angles a, (3, 7 which specify orientation 
of the molecular frame with respect to the laboratory one. Let us further introduce the 
conditional probability density function, 

P (n \n,t), (16) 



which is the probability density that the molecule has orientation f2 at time t, provided it 
had orientation Qq at t = 0. By definition, the quantity f|T6l) obeys the initial condition 
p(O |fi, 0) = 5(0 — O ). We also assume that the molecule can be subjected to an external 
(anisotropic) potential U(fl), so that the corresponding equilibrium Boltzmann distribution 
reads 

pu{0) = Zuexp{-U(n)/(k B T)}, (17) 

ks being the Boltzmann constant, T being the temperature, and Zy being the partition 
function. Evidently, p(O |f2,t — > oo) = p c/ (f2). 

We are in a position now to define the OCF— iM<&L 

G j l, mm ,(t) = (&um)DlUm)) , (is) 

Dj m (H) being the Wigner D-functions.— If the conditional probability density function (fl6|) 
is known, then the OCF fflHl) can be evaluated as follows: 

G% mm >{t) = J dndn Dijn)Dfc,(n )pu(n )p(n \n,t). (19) 

The OCFs ( !T8|) have been explicitly computed, for example, within the rotational diffusion 
model.— ^ 

If the OCFs ( fl8|) are known, one can evaluate any polarization response of interest. 
Indeed, the CF of any (for simplicity, real) orient at ion- dependent quantity B(Q) at the 
time moments and t can immediately be expressed through OCFs ( fl8l) : 

oo j j' 

<s(n(o))B(n(f))> = £ £ (20) 

j,j'=0 l,m=-j l',m'=-j' 

Here 

B L - J dn D jjn)B(n) (21) 

are the expansion coefficients of the quantity B(Q) over the D-functions. 

There exists an important particular case, in which CF (1201) simplifies greatly. Namely, 
let us assume that there are no external fields (U(n) = 0), so that rotational phase space is 
isotropic. Then one can show (see, e.g., Ref.— for more details), that the fundamental OCF 
( JT8l) becomes 

G^mm'W = 2 • + ^ 3i' S U'Grnf m (t) f (22) 
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8jj> being the Kronecker symbol. Therefore, the CF (I20I) is now evaluated as 

oo j 

(B(fl(0))S(n(t))> = E 4^<&m(*), (23) 

j=0 m,m'=—j 

I 3 ■ 

^mm' = 2j _j_ i E ^km^im 1 (24) 
fe=-j 

(compare with^ 1 ^ 1 ^). The quantities 

^w(t) ss £ (DL k m))DZ k (m)) = (D mm '(^(t))) • (25) 

k=-j 

are the standard OCFs of the rank j (fl^(t) being the angle of relative reorientation). 
They can either be simulated on a computer— or evaluated within the models of molecular 
reorientation available in the literature.- 1 '^ 

Eq. (|23l) can further be simplified in the following important particular case. Let us 
assume that the quantity of interest, B(Q), is specified by the orientation of a unit vector 
(for example, a dipole moment) d, which is "rigidly attached" to the molecule, i.e., B(Tt) — > 
B(d). The orientation of the unit vector d is uniquely determined by the the spherical 
angles a, (3 and is specified, therefore, by the reduced D-function Dj (a, (3, 0) which, apart 
from the numerical factor, coincides with the corresponding spherical harmonics,— 

Y jl (a,(3) = ^^Di (a,(3,0) (26) 

On writing Eq. (|26p . we have tacitly assumed that the unit vector d is pointed along the 
z-axis of the molecular reference frame. In this case, the terms with m = ml = survive 
only in Eq.( l23i) . so that 

oo 

(2J(d(0))B(d(f))> = £ Al Gl (t). (27) 

3=0 

Thus, from the formal point of view, the CF ([27]) is determined by the linear combination 
of the OCFs G J Q0 (t) with different ranks j. Of course, the relative significance of the contri- 
butions from different OCFs is determined by the weighting coefficients t4q , which depend 
on a particular form of the quantity under study, B(d). 



7 



IV. SINGLE MOLECULE SIGNAL THROUGH ORIENTATION AL CORRELA- 
TION FUNCTIONS 



The CF, which is normally extracted from the SM polarization signal, is defined via the 
expression2^^ii2^i2i 

c(t) = (x(n(p))x(n(t))), (28) 

in which X(Q) is the SM dichroism f|T2|) and its time-dependence is induced by molecular 
rotation.— To apply a general formalism outlined in the previous section for the evaluation 
of Eq. (|28|) . one has to express all the scalar products in Eq. (|12|) in terms of the Wigner 
D-functions. This is easily achieved by the formula 



M 2 = Ml + 2 J2 ^(0,-/? Sl -a s )4(O)^K,ft,0) 

V l,m=-2 , 



(29) 



and similar expressions for the scalar products (e^/i) 2 . Emphasize that the spherical angles 
a s ,j3 s and a M , (3^ in Eq. fT2T?l) are time-independent and specify orientations of the vectors 
s and fi, correspondingly, in the laboratory and molecular reference frames. The proper 
description of molecular reorientation is accounted for by the (time-dependent) Euler angles 

n. 

To simplify the subsequent calculations, we can proceed as has been described in the 
previous section and choose the molecular reference frame in such a way that the emission 
dipole moment [i is directed along the z-axis of this frame. This is tantamount to putting 
a,j = ^ = in Eq. (l29l) . Then D J m0 (0,0, 0) = 5 m o (the latter being the Kronecker delta), 
and Eq. (1291) becomes independent of the angle 7 of rotation about the molecular z-axis: 



M 2 = ~ ^1 + 2 ^(0, ~Ps, -a s )Df (a, (3, 0) 



(30) 



Upon the insertion of this formula (and similar expressions for (ej/x) 2 ) into Eq. (121]) one 

realizes that dichroism X (Eq. fll2l) ) becomes a function of two Euler angles a and /3, X(a, (3). 

Thus, Eqs. (T21~j) - (T25|) are applicable to this case. By using the interconnection ff26l) between 

the spherical harmonics and D-functions, one finally gets 

00 j 
C(t) =J2 A ooGUt), Ai = l^ol 2 > (31) 

j=0 k=-j 
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/ 1 f n f 

Xi =^—J sm(3d(3 J daY jh (a,P,0)X(a,P). (32) 

These are exactly the formulas which have been derived in Refs.— ^ Emphasize that the 
domain of validity of Eqs. (l3Tj) . (132]) is not confined to linear rotors and spherical tops. 
They are valid for a general asymmetric top molecule and will be employed in the remainder 
of the present paper. 

A word of caution is however in order. First, one should keep in mind that Eqs. (EH), 
(I3"2"|) are applicable for molecular rotation in an isotropic space. If there exist external 
potentials, one should use a more general formula (|2"Uj) . Second, Eqs. (I3~lj) . (1521) are valid if 
and only if the molecular z-axis coincides with the direction of the emission dipole moment 
/i. This direction, however, may not coincide with the molecular symmetry axis. This means 
that such a choice of the molecular frame may not accommodate the molecular symmetry 
properly. In order to do that and to evaluate G J 00 (t) efficiently, we can switch from the 
initial frame to another "convenient" molecular frame via the corresponding Wigner matrix 
D J mk (—A.), evaluate the OCF G J kn (t) in the reference frame which fully accounts for molecular 
symmetry and return back to the original molecular frame: 

GUt) = £ Di k (-A)G{ n (t)Di (A). (33) 

k,n=—j 

Emphasize that the angles A which specify the relative orientation of the molecular frames 
introduced above are fixed and known for any particular molecule. To illustrate the use 
of Eq.( l33l . let as consider a perpendicular transition of a symmetric top molecule. Then 



G L = 6 L G L and 

j 

GUt) = ]T DUO, -vr/2, 0)GUt)D j n0 (0, vr/2, 0). (34) 

n=-j 

For the small angle rotational diffusion, for example, 

G J nn = exp{-j(j + l)D\\t - n 2 (D ± - D\\)t}, (35) 

D\\ and D± being the corresponding diffusion coefficients.— Asymmetric top OCFs G J 00 (t) 
within the small angle rotational diffusion equation can be computed, e.g., by the method 
described in.— 

For the sake of the further comparison, we also present the standard formulas for the 
intensity of the polarized emission in ensemble measurements.- 1 ^ Incorporating the cross- 
sections ([3]) and into Eq.(T2]) and applying Eqs. ff23]) - ff25]) . one obtains the standard result 



that the averaged emission intensity is uniquely determined by the second-rank OCF of the 
dipole moments involved, 

(h) ~ 1 + ^P 2 (e tt eO (P 2 ^aKt)) • (36) 

Here Pj(x) are the j-rank Legendre polynomial and 

i 

{PjMt)) = E C^-a^-^JC'W^oK.^-O), (37) 

m,m'=—j 

a na> flva an d a M' ^ being the spherical angles which specify the orientation of the absorption 
(Ha) and emission (/x) dipole moments in the molecular reference frame. 

V. RANK DEPENDENCE OF THE SINGLE MOLECULE SIGNAL 

Once the emission dipole moment \x is pointed along the z-axis of the molecular reference 
frame, then Eqs. (13"TT) and (1521) determine the SM CF. The coefficients A J 00 depend upon 
the OCF rank j, the NA angle x-, an d the relative orientation <p of the polarizers ei and e 2 
{if = arccos(e!e 2 )). Hereafter, we denote these coefficients as <p). The standard choice 
of the orthogonal detection scheme corresponds to A^(x, 7r/2). The term A°(x, ip) = 0, since 
the isotropic component does not contribute to dichroism ( fl2l) . Furthermore, the symmetry 
of the D- functions dictates that only the coefficients with even j contribute into Eq. (l3Tl) . It 
can thus be recast into the form 

oo 

C(t)=$> 2CT ( X ,y,)G 2CT (t). (38) 

CT=1 

(we have dropped the 00-subscripts for the clarity of notation, and the CF is assumed to be 
normalized to unity, that is C(0) = 1). Formally speaking, the SM CF is expressed by the 
linear combination of all even-rank OCFs.— Their relative contributions are determined by 
the weighting coefficients A 2a (x, ip), which decrease rapidly with the OCF rank j, so that a 
few of them contribute significantly into Eq. (l38l . 

Let us consider a special case of the orthogonal registration scheme (<p = vr/2) first. If the 
NA is negligibly small (x = 0), one gets the set of coefficients A 2r7 (0, 7r/2), which have been 
calculated and analysed in.- 2 ^ I n this case, A 2 (0,vr/2) = 0.835, A 4 (0,n/2) = 0.100, etc., so 
that the contribution due to the second-rank OCF yields more than 83%. In the opposite 
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case of a high NA (x = 7r /2), as is clear from Eq. (fl5|) . the second-rank OCF contributes 
exclusively into CF fl38|) . and C(t) ~ (P 2 (/i/i(t)). Thus, apart from the numerical factor, 
the SM CF coincides with the ensemble averaged anisotropy. Since the parameters p and 
q in Eq.flSj) increase monotonously with \ ( see Fig- 2), the significance of all coefficients 
A 2a (x, 7r/2) with a > 1 decreases with x- As has been shown in,— for example, the CF C(t) 
corresponding to NA = 0.6, n = 1.4 coincides, practically, with G 2 (t). 

The analysis which has been carried out in Ref.-^ and the above considerations mean that, 
under some typical experimental conditions (ip = n/2 and x being close to 7r/2), the SM CF 
C(t) is indistinguishable from the second-rank OCF G 2 (t). If one considers a complicated 
dependence of C(t) on OCFs of various ranks as a nuisance, which obscures interpretation 
of the measured signal, then the use of the orthogonal detection scheme in conjunction with 
high-NA optics allows one to "measure" the standard second-rank OCF. On the other hand, 
Eq. (l38l) hints at a unique opportunity to learn about higher-rank OCFs experimentally, even 
with the use of high-NA optics. We suggest that the non-orthogonal registration scheme, 
when ei and e 2 are not perpendicular to each other, makes this possible. 

To clarify the situation, let us consider the signal intensity ( 1361) which is measured within 
the ensemble- averaged spectroscopy. The intensity consists of the sum of the isotropic 
component, which contains no information about molecular rotation, and the anisotropic 
component, (P 2 (/x a /z(t)). The effect of the detection scheme is exclusively defined by the 
numerical factor P2(e a ei), which determines the relative weights of the two components. 
Therefore, in order to extract OCF (P 2 (/i a /i(t)) out of the signal, it is sufficient to perform 
measurements of the emission intensities at any two different polarizations ei and e 2 . A 
common practice is to use the magic angle conditions (P 2 (e a ej) = 0, e a ej = 1/V3), as well 
as parallel (P 2 (e a e,) = 1, e a ej = 1) and perpendicular (P 2 (e a ej) = —1/2, e a ej = 0) detection 
schemes. 

The situation with the SM polarization spectroscopy is very different. The SM CF f )38|) 
has weighted contributions from virtually all even-ranked OCFs, and the weights themselves, 
A 2a (x, <f), are clearly detection scheme dependent. The coefficients A 2cr (x, <p) for several few 
first a are plotted in Fig. 3 for % = (small NA) and for % = ir/2 (high NA). The 
parameters p and q in Eq. (|12p increase monotonously with x (see Fig. 2). One thus expects 
that the coefficients A 2cr (x, <p) transform gradually from those depicted in Fig. 3 in the 
upper panel to those depicted in the lower panel, following the increase of x- 
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Evidently, A 2 (x, ^-components dominate the signal for any polarization scheme, so that 
the contribution due to the second rank OCFs is the most significant. Both A 2 (0,ip) and 
A 2 (n/2, ip) reach their maxima at ip = ir/2: A 2 (0, <p> = ir/2) = 0.835 and A 2 (ir/2, ip = tt/2) = 
1. This means that the standard orthogonal scheme reflects predominantly the decay of the 
second rank OCF, since 84% of the low-NA signal and 100% of the high-NA signal is 
determined by G 2 (t). If the angle ip between the polarizations ei and e2 decreases, then the 
contribution due to A 2 (0,ip) also dominates, but the higher rank OCFs start to contribute 
more and more significantly. For a \ow-NA signal, for example, the forth-rank contribution 
achieves its maximum of 23% at tp = it/ 8 while the second order contribution remains as 
high as 46%. The shapes of the curves A 2a (0,ip) and A 2r7 (7i/2, ip) are evidently different. 
Every coefficient A 2cr (0, ip) with a > 1, as a function of (p, exhibits an (asymmetric) bell-like 
shape with a single maximum. A 2 (tt/2, ip), as well as A 2 (0, ip), increase monotonously with 
ip, while all A 2a (ix/2, ip) with a > 1 decrease rapidly. An overall tendency can be described 
as follows: the closer are ei and e2 to each other, the more A 2a (x, <p) coefficients contribute 
to the signal. This is clarified by Figs. 4, in which we present the completeness parameters 



for few first £. As to the low-NA detection, the first few a are necessary to faithfully 
reproduce the signal for ip close to tt/2, while much more terms are necessary for small ip. 
The same is also true for the high-NA detection, but the convergence is much more rapid. 
In general, the high-NA detection is more second-rank OCF dominated than the low-NA 
one. We emphasize however that the high-NA CF C(t) coincides with G 2 (t) in the case of 
the orthogonal detection scheme only (<p = tt/2). Otherwise, it also depends on OCFs of 
different even ranks. 

We conclude the present section with the following comment. On writing the starting 
Eq.([T]) for the SM dichroism, we have tacitly assumed that the emission intensities I\ and 
I2 scale identically. To take into account a possible imbalance rj of the channels, we can 
redefine the SM dichroism as (ii — ^7-^2) /(^l +V^)-~ We can then repeat the above analysis, 
taking into account that the coefficients A 2a (x, ip) acquire an additional ^-dependence. In 
that case, for example, the isotropic contribution is nonzero, A°(x, f) ~ 1 — 77 7^ 0. If, 
furthermore, there are several emission dipole moments,™ the above approach can also be 
generalized straightforwardly. 




(39) 



cr=l 
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VI. EXTRACTION OF THE HIGH-RANK ORIENTATION AL CORRELATION 
FUNCTIONS FROM THE SINGLE MOLECULE SIGNAL 

As has been demonstrated in Sec. 5, the non-orthogonal detection schemes contain, po- 
tentially, more information on the high-rank OCFs than the standard orthogonal schemes. 
The question thus arises if it is possible to extract G 2cr (t) with a > 1 from the SM CF 
C(t). To clarify the situation, we compute the signal within the simplest model of molecular 
reorientation, the spherical top small angle diffusion. Within this model,— the OCFs are 
described by the exponential formula 

Gi (t) = GUt) = exp{-j(j + l)D l{ t}. (40) 

We have chosen this model since it is most frequently applied for the interpretation of SM 
CFs. 

Let us consider the small- NA detection first (Fig. 5). C(t) calculated within the orthog- 
onal detection scheme (ip = vr/2) and within the scheme with ip = tt/8 (this angle provides 
the maximum for the G 4 (t)-contribution, see Fig. 3, upper panel) are seen to be markedly 
different, both mutually and from G 2 (t). On the other hand, the SM CF with p = n/2 
is rather close to G 2 (t), and the slopes of both SM CFs, as well as the slope of G 2 (t), are 
almost the same. This is totally understandable, since (i) 84% of the p = ir/2 CF and 46% 
of its ip = tt/8 counterpart are determined by G 2 (t) and (ii) the higher-rank OCFs decay 
much more rapidly than those with j = 2. 

Theoretically speaking, the procedure of the extraction of the high-rank OCFs from C(t) 
is straightforward. One can perform several (say, N) measurements at N different detection 
angles ip, truncate the number of summations in Eq. (138|) by N, and solve the corresponding 
system of linear equations for G 2a (t), a = 1 4- N. Such a procedure, however, can hardly 
be feasible in practice. There exists a cruder, but much more robust procedure which is 
exemplified by Fig. 5. We can calculate the quantity 

C(t) = C(t, <p = tt/8) - * = n/2) - (41) 

Evidently, C(t) represents a weighted sum of the OCFs of the rank j = 4 and higher. Since 
the higher-rank contributions into C(t) decrease rapidly with a (Fig. 3) one expects that 
C(t) is determined, predominantly, by G A {t). This qualitative expectation is corroborated 
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by Fig. 5. Evidently, C(t) and G 4 (t) do not coincide but, as has been explained above, their 
slopes are virtually the same. 

The procedure described above can also be applied to the high-iVA detection, see Fig. 6. 
Furthermore, the situation is much more fortunate in this case, since the orthogonal CF is 
now solely determined by the second rank OCF G 2 {t) (A 2 (n/2, ir/2) = 1, see also Fig. 3). It 
is interesting that CFs C(t,tp = 7r/4) and C(t, (p = 7r/2), which are presented in the Fig. 6, 
are rather close to each other, since both of them are predominantly determined by G 2 (t). 
However, the CF 

C{t) = C{t, tp = tt/4) - A 2 {n/2, n/4)C{t, p = tt/2) (42) 

is almost indistinguishable from G 4 (t). 

The above results demonstrate that the proposed (or similar) procedure of the extraction 
of high-order OCFs from the SM CF C(t) is useful and robust, both for high- and low NA 
objectives. Once the higher order OCFs are available, one can get valuable information 
on the dynamics of the SM reorientation. For example, let us suppose that orientational 
relaxation proceeds exponentially. Then, if the second-rank OCF is available only, one can 
extract an effective time of the OCF decay, but cannot learn anything about rotational 
dynamics. If both G 2 (t) and G A (t) are known, one then can calculate the ratio of their 
relaxation times. This quantity, being very model-specific, allows one to discriminate be- 
tween different reorientation mechanisms. For the small-angle diffusion model, the ratio 
equals to 10/3 (see Eq. (l4"0l ). For the jump diffusion model, it is close to i.i&2L2& More 
sophisticated approaches to the orientational relaxatioi^ 1 ^^ 1 ^ 1 ^ 1 ^ 1 ^ also predict that the 
comparison of the behavior of OCFs of different ranks makes it possible to identify the 
underling mechanisms of orientational relaxation. 

Before concluding the present section, we wish to discuss possible causes of the deviation 
of C(t) from the exponential form (see also Refs.— ^). There can be two fundamentally 
different groups of reasons. First, C(t) is described by the weighted sum of the even-rank 
OCFs, Eq. (l38l) . Despite the second-rank OCF contributes predominately into C(t), the 
contributions due to the higher-rank OCFs cannot be neglected, in general. It is only for a 
high-NA objective and close-to-orthogonal detection scheme that these higher order contri- 
butions are vanishingly small and C(t) ~ G 2 (t). Second, the OCFs (including the second- 
rank OCF) can be non-exponential due to a variety of reasons, (i). Molecular rotation is 
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not necessarily diffusive. While the jump diffusion model predicts exponentially decaying 
OCFs^^ 1 ^ (although their rank-dependence differs from that given by Eq. fl35|) ). the re- 
stricted diffusion model,— ^ the diffusion-equation-with-memory models and other memory 
function approaches^^ 1 ^ predict the OCFs of a spherical molecule to be described by the 
sum of several (in general, complex) exponentials. More sophisticated approaches deliver, of 
course, more complex OCFs. On the other hand, the so-called inertial effects, which induce 
highly non-exponential behavior of OCFs,->^>2£ are irrelevant for the SM spectroscopy since 
they manifest themselves on a time scale which is much faster than the time resolution of 
typical SM experiments, (ii). A deviation of the molecular shape from spherical complicates 
molecular reorientation even in the small angle diffusion limit. The second-rank OCF Gqo(^) 
of an asymmetric top, for example, is described by the sum of two exponentials.— 1 ^ (iii). Of 
relevance is the direction of the emission dipole moment /i in the molecular reference frame 
(see Eqs. (1331) and For example, if /i is parallel to the axis of the linear or symmetric 

rotor, the corresponding second-rank OCF is mono-exponential. If /i is perpendicular to the 
molecular axis, then the corresponding OCF is two-exponential (see Eq. ([35l) ). (iv). Finally, 
internal rotations can also cause deviations from exponentiality. 46 

VII. CONCLUSION 

We have developed an approach for the calculation of the SM CF C(t), Eq. (T25i) . in 
the general case of asymmetric top molecules. The key dynamic quantities of our analysis 
are the even-rank OCFs G^(t) fl25|) . the weighted sum of which constitutes C(t). The 
OCFs can either be evaluated within any model of molecular reorientation available in the 
literature^^ Q r simulated on a computer.— We have demonstrated that the use of non- 
orthogonal schemes for the detection of SM polarization responses allows one to manipulate 
the weighting coefficients in the expansion of C(t) on OCFs. Thus a valuable information 
about the OCFs of the rank higher than second can be extracted from the SM CF C(t). 
Until recently, such an information was available only through computer simulations and/or 
model calculations. Neither the corresponding information is accessible within the ensemble- 
averaged spectroscopy, in which one "measures" exclusively the second-rank OCF (l36l) . 
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FIGURE CAPTIONS: 



FIGURE 1: Schematic representation of the SM polarization experiment. 
FIGURE 2: The iVA-angle dependence of the coefficients p (dashed line) and q (full 
line) . 

FIGURE 3: The expansion coefficients A 2a (0,p) (low-NA objective, upper panel) and 
A 2cr (7r/2, ip) (high-NA objective, lower panel) as functions of the angle p> between the polar- 
izers. From top to bottom, the curves correspond to a = 14-8 (upper panel) and a — 1 4- 6 
(lower panel). 

FIGURE 4: The completeness parameters B 2a (0, (p) (low-iVA objective, upper panel) 
and B 2r7 (n/2,(p) (high-NA objective, lower panel) as functions of the angle p> between the 
polarizers. From top to bottom, the curves correspond to o = 1 4- 8 (upper panel) and 
a = lv4 (lower panel) . 

FIGURE 5: Extraction of the higher-order OCFs from the single molecule CF C(t) (|38|) 
in the case of a low-NA objective (x = 0). The upper solid curve depicts the the second-rank 
OCF G 2 (t) and the lower solid curve depicts the forth-rank OCF G A (t). The upper dashed 
curve shows the single molecule CF C (t) calculated for p = tt/8 and the lower dashed curve 
shows the single molecule CF C(t) calculated for p = n/2. The dotted curve depicts the 
forth-rank OCF, which is approximately extracted from the two single molecule CFs (see 
text for details). All the CFs are calculated within the small angle rotational diffusion model 
for a spherical top, D\\ =0.1; t and D^ 1 are given in arbitrary dimensionless units. 

FIGURE 6: Extraction of the higher-order OCFs from the single molecule CF C(t) 
(13"8"|) in the case of a high-NA objective (x — 7r /2)- The upper curves depict the second- 
rank OCF G 2 (t) (solid line) and the single molecule CF C(t) calculated for p = n/4 (dashed 
line). The lower curves show the exact forth-rank OCF G 4 (t) (solid line) and its counterpart 
extracted from CF C(t) (dashed line, see text for details). All the CFs are calculated within 
the small angle rotational diffusion model for a spherical top, D\\ = 0.1; t and D7 are given 
in arbitrary dimensionless units. 
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